Role of Ponded Turbidity Currents in Reservoir Trap Efficiency

The capacity to store water in a reservoir declines as it traps sediment. A river entering a reservoir forms a prograding delta. Coarse sediment (e.g., sand) deposits in the fluvial topset and avalanching foreset of the delta, and is typically trapped with an efficiency near 100%. The trap efficiency of fine sediment (e.g., mud), on the other hand, may be below 100%, because some of this sediment may pass out of the reservoir without settling out. Here, a model of trap efficiency of mud is developed in terms of the mechanics of a turbidity current that plunges on the foreset. The dam causes a sustained turbidity current to reflect and form a muddy pond bounded upstream by a hydraulic jump. If the interface of this muddy pond rises above any vent or overflow point at the dam, the trap efficiency of mud drops below 100%. A model of the coevolution of topset, foreset, and bottomset in a reservoir that captures the dynamics of the internal muddy pond is presented. Numerical implementation, comparison against an experiment, and application to a field-scale case provide the basis for a physical understanding of the processes that determine reservoir trap efficiency.     

Practical Equivalent Continuum Model for Simulation of Jointed Rock Mass Using FLAC3D

A simple practical equivalent continuum numerical model for simulating the behavior of jointed rock mass has been extended to three-dimensional using FLAC3D. This model estimates the properties of jointed rock mass from the properties of intact rock and a joint factor (Jf), which is the integration of the properties of joints to take care of the effects of frequency, orientation, and strength of joint. A new FISH function has been written in FLAC3D specifically for modeling jointed rocks using the Duncan and Chang hyperbolic model. This model has been validated first with simple element tests at different confining pressures for different rocks with different joint configurations. Explicit modeling of the joints has also been done in element tests and results obtained compare well with the results of equivalent continuum model and also with experimental results. Further, this has been applied for a case study of a large underground power house cavern in the Himalayas. The analysis has been done under various stages of excavation, assigning a null model available in FLAC3D for simulating the excavation.     

Applications of the discrete element method in mechanical engineering

Compared to other fields of engineering, in mechanical engineering, the Discrete Element Method (DEM) is not yet a well known method. Nevertheless, there is a variety of simulation problems where the method has obvious advantages due to its meshless nature. For problems where several free bodies can collide and break after having been largely deformed, the DEM is the method of choice. Neighborhood search and collision detection between bodies as well as the separation of large solids into smaller particles are naturally incorporated in the method. The main DEM algorithm consists of a relatively simple loop that basically contains the three substeps contact detection, force computation and integration. However, there exists a large variety of different algorithms to choose the substeps to compose the optimal method for a given problem. In this contribution, we describe the dynamics of particle systems together with appropriate numerical integration schemes and give an overview over different types of particle interactions that can be composed to adapt the method to fit to a given simulation problem. Surface triangulations are used to model complicated, non-convex bodies in contact with particle systems. The capabilities of the method are finally demonstrated by means of application examples. Commemorative Contribution.     

Numerical solution for shape optimization of an impactor penetrating into a semi-infinite target

The shape of the impactor with the maximum depth of penetration (DOP) for a given impact velocity is found using a numerical procedure for solving a corresponding non-classical variational problem. It is shown that the optimum shape in a general case is close to a blunt cone. The variation of the optimal shape of the impactor and the dependence of the DOP vs. the initial (impact) velocity and friction coefficient is studied. The analysis is performed also for optimal conical impactors.     

粉尘薄膜爆轰

与普通爆轰不同，粉尘薄膜爆轰波前氧化剂与燃料处于分离状态．由于激波作用下壁面燃料的注入才于壁面附近形成燃烧的粉尘云，燃烧释放的热使激波自持形成粉尘薄膜爆轰。本文在长2、8m，横截面20mm×20mm的方型水平激波管中对此类爆轰进行了实验研究．并根据作者建立的模型对其二维结构进行了数值模拟。     

熔窑玻璃液流数值模拟的新发展

熔窑玻璃液流数值模拟的新发展赵国昌（清华大学核能技术设计研究院，１０００８４）ＮｅｗＤｅｖｅｌｏｐｍｅｎｔｏｆｔｈｅＮｕｍｅｒｉｃａｌＳｉｍｕｌａｔｉｏｎｆｏｒＦｌｏｗａｎｄＨｅａｔＴｒａｎｓｆｅｒｏｆＧｌａｓｓＭｅｌｔｉｎＴａｎｋＦｕｒｎａｃｅ￥Ｚ...     

NUMERICAL SIMULATIONS OF WHOLE FIELD BEHAVIOR FOR ROUND BUOYANT JET

ＮＵＭＥＲＩＣＡＬＳＩＭＵＬＡＴＩＯＮＳＯＦＷＨＯＬＥＦＩＥＬＤＢＥＨＡＶＩＯＲＦＯＲＲＯＵＮＤＢＵＯＹＡＮＴＪＥＴＬｉＷｅｉ，ＨｕａｉＷｅｎ－ｘｉｎ（ＷｕｈａｎＵｎｉｖｅｒｓｉｔｙｏｆＨｙｄｒａｕｌｉｃａｎｄＥｌｅｃｔｒｉｃＥｎｇｉｎｅｅｒｉｎｇ，...     

油砂地层中子-伽玛能谱Monte-Carlo模拟

本工作应用自行研制开发的中子 伽玛能谱测井的中子 伽玛射线随空间、能量、时间分布的蒙特卡罗模拟软件包计算了裸眼井、井眼里无下井仪器、油饱和砂岩地层、井眼注油条件下沿井轴和井壁的中子 伽玛射线随能量、时间分布及沿井轴和沿井壁的中子 伽玛射线能谱 ,研究了源距对中子 伽玛射线随能量、时间分布以及井轴与井壁中子 伽玛射线能谱的影响。